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As  a  preliminary  we  examine  a  two -ray  condition  as  shown  in  Fi.'^ure  4  . 
One  wavefront  is  assumed  to  arrive  ale ng  the  x-axis,  another  at  an  angle  0, 
The  wavefronts  are  sinusoidal  in  time  and  the  sum  at  points  along  x  is 

s(t,x)  *=  Aj^cos(ajt  -  kx)  +  A^cos^wt  ~  (f  -  kx  cos0)  (1) 


FIGURE  4,  RAY  ARRIVALS 


Tlie  case  treated  here  is  analogous  to  oac  encountered  in  FM  systems  with 
slnu*'oidal  interference  to  a  desired  carrier. 
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and  represent  the  magnitude  of  the  two  wavefronts,  (J  is  the  phase  dif¬ 
ference  between  them  at  x  •=  0,  and  k  =  27r/X  where  X  is  the  wavelength.  A 
sensor  at  point  x  responding  equally  and  linearly  to  both  wavefronts  would 
see  s(t,x).  Writing  (1)  in  envelope-angle  form. 


where 


s(t,x)  -  B(x)cos(u)t  +  H'(x)) 


B(x) 


2A1A2 


cos  [(J)  -  kx(l  -  COS0)] 


i|;(x)  =  tan 


-1 


A^sin  kx  +  A2sin(<)  +  kxcos0) 
Aj^cos  kx  +  A2Cos(4>  +  kxcos8) 


(2) 

(3) 

(4) 


The  phase  obtained  using  (4)  will  be  inodulo-2T7.  It  is  useful  to  deal  with 
the  phase  derivative  d\f*/dx  if  the  inodulo-27r  ambiguity  is  to  be  avoided.  It 
can  be  shown  that 


dx 


k  -  ^(1  -  cos9) 


where 


2 

1  +  a  +  2a  cos[({)  -  kx(l  -  cos0)  ] 


(5) 


a 


A1/A2 


(6) 


When  ^  is  large,  meaning  that  the  important  part  of  the  received  vrave  is 
along  the  x-axis 


=  k  -  cos[(}>  -  kx(l  -  cos9)],  (7) 

It  fluctuates  sinusoidally  around  k  with  the  fluctuation  amplitude  decreasing 
to  zero  as  ^  goes  to  infinity.  For  ^  small 

=  k  COS0  -  ka(l  -  cosO)  cos[(J)  -  kx(l  -  cos0)],  (8) 

ngaln  a  sinusoidal  fluctuation  which  decreases  to  zero  as  ^  approaches  zero. 
Fol  intermediate  values  of  the  fluctuation  of  di»/dx  Is  as  shown  in  Figure  5, 
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FIGURE  5 .  PHASE  DERIVATIVE  VARIATION  WITH 
DISTANCE  IN  THE  TWO-RAY  CASE. 

The  variation  is  periodic  with  period 


X  = 


1  -  cos© 

and  has  peak  excursions  above  and  below  k  given  by 


(9) 


g  «  k(l  -  cos0)/(a  -  1)  (10) 

a  =  -k(l  -  cos0)/(a  +  1)  (11) 


Figure  5  is  drawn  assuming  a  >  1.  The  fluctuation  is  around  k,  which  turns 
out  to  be  the  average  of  di|>/dx.  Note  that  at  a  *=  1^  B  is  positive  and  high 
in  magnitude,  and  the  fluctuation  is  highly  impulsive.  For  a  <  1,  di(;/dx 
fluctuates  around  the  value  k(l-cos6)  rather  than  around  k  and  for  a  =  1  B  ir 
negative  and  high  in  magnitude,  and  the  fluctuation  is  again  impulsive  but 
negative  going.  When  ^  is  close  to  unity  the  phase  as  a  function  of  posi¬ 
tion  (which  is  the  integral  of  dij//dx)  is  as  shown  in  Figure  6, 

In  underwater  applications  the  angle  0  typically  found  in  long  range 
paths  is  less  than  20^.  Assuming  it  to  be  10^  the  period  is 


X 


\ 

1  -  0.98 


50X 


(12) 


> 
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It  Is  interesting  to  note  that  in  analyses  found  in  the  literature  of  the  cor 
relation  distance  of  underwater  acoustic  waves,  numerical  estimates  around 
SOX  are  typically  obtained  (see  for  Instance  II)). 


tl]  P.  VT.  Smith,  Jr.,  ’’Spatial  Coherence  in  Multipath  oi  Hultimodel  Channels,*’ 
Journal  Acoustic  Society  of  Amerlra,  Vol.  60,  No.  2,  Aug.  1976,  pp.  305-310 


We  now  turn  to  the  more  general  case  of  n  wavefronts  arriving  at 
angles  0^,  j  **  l,2,,.>.n.  The  received  sum  at  a  position  x  is 


s(t,x)  »  I  A.cos(a)t  -  -  kx  cos0.) 

j  j  j 


(13) 


The  j  «=  l,2,...n  are  random  phase  angles  of  each  of  the  wavefronts  on 
arrival  at  the  point  x  =  0.  It  is  convenient  to  write  this  in  the  form 


s(t,x)  *  B(x)  cos[a)t +il;(x)  1 


jmt 


where 


■■  Re  z(x)e 


"  j  ($  +  kx  COS0  ) 

:(x)  =  B(x)eJ'*'^''^  =  E  A  e  ^  J 

j=l  J 


(14) 


(15) 


iKx)  is  the  phase  angle  we  will  study  ar.d  as  was  done  before  we  find  the 
phase  derivative 


di|)  _  ^  ,  1 


From  (15)  we  have 


,  n  j(<|i.  +  kx  COS0  ) 

-J-  *=  j  I  A.  k  cosO,  e  ^  ^ 

Ox  ■'j.i  J  j 


SO  that  (16)  becomes 


dx 


Im  j 


j  (4> .  +  cosO  ) 
lAj  k  COS0J  e  ^  ^ 


1 


EA.  e 

J  ^ 


j  ((J)j  +  kx  cosOj) 


E  E  A  A,  k  cosO.  c 
i  i  ^ ’ 


j[<'.  “  kx(cosO^  -  cosO^)] 


E  A.  c 


cosO . )  I  ? 


E  E  ^  cosO^  +  kx(cosO^.  -  cosOj)! 

E  E  +  kx(c()sO^  -  ro  ;0j)l 


(16) 


(17) 


(18) 


OVR  N^. 


1 


Where  this  approximation  is  permis- 


As  a  rule  9^  «  cosG^  == 

sible  ve  can  write 


dx 


k 


k 

2 


Z  Z  0^^ 

Z  Z  a^j(x) 


where 


aij(x)  =  A^Aj  cos[(f^  ~  ^  kx(cos0^  -  cosO^)] 


To  retrieve  i{»  we  must  integrate  di^i/dx;  l.e.  we  form 


(19) 


d,p(x  ) 

'  dx~~  °  ~ 

'o 


The  integration  will  generate  the  phase  difference  between  the  phase  at  x 
and  the  phase  at  the  origin  of  integration.  From  (19)  we  see  that  one  term 
on  integration  will  te  kx,  the  linear  phase  variation  associated  with  the 
normal  phase  vs.  position  function  of  a  plane  wave  along  the  direction  of 
travel  of  the  wave.  In  beam  forming  with  an  array  of  sensors  along  x  one 
will  subtract  the  phase  progression  kx  if  the  axis  of  the  beam  is  to  be 
colinear  with  the  x  axis.  In  this  case  the  remaining  phase  difference  be^ 
tween  a  point  x  and  the  origin  is 


-k 

2 


Z  Z  0,^  a 

LJ 

Z  Z  a 
o  i  j 


Ij(Xi) 


Ij(Xi) 


dx. 


(21) 


2 

If  the  approximation  cosOj^  «  1  -  9^  /2  is  not  used*  the  renviinlng  phase 
after  correcting  for  kx  Is  given  by  subtracting  k  from  (18)  and  integrat¬ 
ing  over  X. 

Numerical  evaluations  of  the  remaining  phase  difference  have  been  made 
for  a  number  of  cases. One  parti cxtlar  case  is  shown  in  Figure  7  determined 
assuming  21  equal  amplitude  rays  arriving  at  2°  Intervals  from  0  =  -20^^  to 


A 

Prograiflinlng  of  this  computation  and  the  one  described  later  giving  array 
pattern,  wat*  done  by  On  dtian  lio. 
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2nrx 


FIGURE  7.  PHASE  DIFFERENCE  V3  POSITION 


I 


+20*  with  respect  to  the  horizontal  (see  Figure  4),  each  with  a  phase  angle 
^  randomly  selected  in  the  interval  (0,2Tr).  Because  the  average  wavelength 
of  the  various  rays  as  seen  along  x  is  less  than  X,  the  wavelength  along  the 
direction  of  travel  of  the  ray,  there  is  a  linearly-tending  phase  accumula¬ 
tion  with  distance  as  seen  in  Figure  7-  On  top  of  this  accumulation  there 
is  a  random  variation.  The  fluctuation  around  a  straight  line  approximation 
to  the  phase  difference  ranges  around  +3  radians.  Thus  even  if  the  phase 
were  corrected  to  account  for  the  slope  of  the  straight  line  approximation, 
a  +  3  radian  random  error  would  still  be  encountered. 

Figure  7  was  obtained  with  one  randomly  selected  set  of  ray  arrival 
phase  angles.  Additional  examples  will  be  ultimately  computed  for  different 
sets  of  arrival  phase  angles  to  provide  data  suitable  for  obtaining  statisti¬ 
cal  averages.  Other  cases,  including  different  intervals  of  arrival  angle, 
different  ray  amplitudes,  and  different  numbers  of  arriving  arrays  will  also 
be  treated. 

Having  a  sample  function  of  phase  vs.  position,  a  logical  next  step  is 
to  determine  the  gain  and  pattern  of  the  random  planar  floating  array  when 
it  is  focused  in  some  azimuthal  direction  using  conventional  beamforming, and 
when  the  source  signal  is  propagating  toward  the  array  through  the  multipath 
medium.  As  a  first  step  a  program  was  developed  for  selecting  element  posi¬ 
tions  over  a  circular  area  assuming  a  uniform  distribution  of  clement  posi¬ 
tions. 

If  the  array  is  assumed  confined  to  a  circle  of  radius  P  with  uniform 
distribution  over  the  circle,  the  density  function  in  the  joint  random  vari¬ 
ables  X,  Y,  is 


!  \ 


Px  y(x,y)  =  — 2 
'  IT  p 


2^2  2 
X  +  y  <  P 


elsewhere. 


Transforming  to  polar  coordinates,  (R,C),  we  have 


■>R,  •  °  1 ' 1  P 

0<<t><2TT 
=•0  ,  elsewhere 
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The  marginal  densities  in  R 

I  \  2r 
Pj^(r)  -  2 

P 

*=  0 

KC  0 


and  ^  arc 

0  ^  r  ^  p 

,  elsewhere 

,  0  <<{.<_  27r 

,  elsewhere 


The  random  variables  R  and  4>  are  independent  and  independent  choices  of  these 
variables  arc  made.  Sample  values  of  ^  are  obtained  by  a  conventional  con*' 
puter  program  which  selects  sample  values  uniformly  distributed  in  (0,1)  and 
multiplies  these  by  2*11.  Sample  values  of  R  are  obtained  by  picking  a  number 
Z  uniformly  distributed  in  (0,1)  and  forming 


P2 


1/2 


for  then 


Pj^(r)  = 


dz 

dr 


0  <  r  <  p  . 


Finally,  the  pairs  (r,cf)  so  obtained  are  converted  back  Into  rectangular  co¬ 
ordinates  by 


X  *=  r  cosef 
y  ”  r  sint^) 

Using  clcnent  positions  so  determined  the  array  pattern  was  next  found. 
The  geometry  of  the  problem  is  sho^n  in  Figure  8, 

Assuming  N  elements  distributed  over  the  circle,  cophased  to  form  a  beam 
along  the  y  axis,  the  array  pattern  is  given  by 

-  N  J[k(x  cos't’  +  y  sin(>  -  y  )  -  a(d  )] 

A(»)  -  I  I  J  BW  )c  "  “  n  n  I 

n»l 


where  x  ,  y  ,  <}>,  and  d  are  defined  in  Figure  8,  and  a(d  )  is  a  phase  vs  po~ 
n  n  n  n 

sltlon  function  of  the  form  obtained  earlier  and  shown  In  Figure  7*  B(d  )  is 

n 

the  amplitude’  of  the  accuistlc  field  at  the.  ntli  elomonl  .  Ihls  qinintily  can  be 
obtained  ur.ing  the  earlier  analyaJe,  but  for  our  pin  posts  now  we  will  assume 
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TYPICAL 

array 

pLEMfKTS 


ORCULAI^  ARRAY 
BOUNDARY, 
RADIUS  . 


FIGURi:  8.  ARRAY  GEOJIETRY. 

it  constant  and  set  it  equal  to  unity  for  all  n.  Amplitude  fluctuations  as 
a  rule,  cause  minor  effects  compared  to  phase  riuctuations .  The  phase  sample 
function  of  Figure  7,  called  now  a(x)  is  used  alone  below  to  assess  the  ef¬ 
fect  of  the  multipath  medium.  The  variable  x  in  Figure  7  is  replaced  by  d  , 

n 

with 

d  =  p  -  r  cos(({)  - 

n  n  a 

f  2  ,  2.1/2  .  --1 

*=  p  ~  (x^  y  )  cos  (tan - v) 

n 

Computer  calculations  of  A(4>),  as  described  above,  were  carried  out  for  two 

cases:  (1)  ^  (2)  ct(d  )  as  given  by  Figure  7,  and  the  results 

are  shown  in  Figures  9  and  10.  Caie  1  is  that  cf  propagation  through  a  trans¬ 
parent  (non-multipath)  medium  while  case  2  is  for  the  particular  multipath 
case  resulting  in  the  phase  function  discussed  above,  Kote  that:  the  gains  along 
the  rnain  beam  in  the  two  cases  are  in  the  ratio  of  about  4.4  dB  -  a  substan¬ 
tial  factor;  the  sidelobe  structure  is  different  in  detail  but  not  in  general 
characteristics.  These  results,  it  imust  be  recognized,  are  based  on  one  set 
of  random  arrival  phases  and  on  on^  set  of  random  element  positions;  whether 
they  are  roprecentativu  remains  tr»  be  detcrinini..d .  Averaging  over  many  sets 
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FIGURE  9  .  ARRAY  PATTER2,'  OF  FLA;;AR  ARRAY  IN  TRANSPARENT 

MEDIUM.  E.NDFIRE  BINV-!  FORMED  AT  ARIJRJTH  AI^'GLE,  =  90®. 


FIGURE  10.  A.RRAY  PATTERN  OF  PYOIDOM  iM.ANAR  AR;l.\Y  ASSUi-lING  FaU.TIPATH 

PROPAGATION.  ENDFIR"  RR':!  iOPJIEU  AT  AElMUFiI  /vNGLE  *  »  90® 
'  FOR  TPv/iNSPARi^NT  MEDILDI. 


of  arrival  phases  and  element  positions,  as  well  as  carryinp,  out  additional 
computations  with  other  system  parameters,  remain  to  be  done. 

Fred  Haber 
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